Evaluation of 16S rRNA Gene Primer Pairs for Monitoring Microbial Community Structures Showed High Reproducibility within and Low Comparability between Datasets Generated with Multiple Archaeal and Bacterial Primer Pairs

The application of next-generation sequencing technology in microbial community analysis increased our knowledge and understanding of the complexity and diversity of a variety of ecosystems. In contrast to Bacteria, the archaeal domain was often not particularly addressed in the analysis of microbial communities. Consequently, established primers specifically amplifying the archaeal 16S ribosomal gene region are scarce compared to the variety of primers targeting bacterial sequences. In this study, we aimed to validate archaeal primers suitable for high throughput next generation sequencing. Three archaeal 16S primer pairs as well as two bacterial and one general microbial 16S primer pairs were comprehensively tested by in-silico evaluation and performing an experimental analysis of a complex microbial community of a biogas reactor. The results obtained clearly demonstrate that comparability of community profiles established using different primer pairs is difficult. 16S rRNA gene data derived from a shotgun metagenome of the same reactor sample added an additional perspective on the community structure. Furthermore, in-silico evaluation of primers, especially those for amplification of archaeal 16S rRNA gene regions, does not necessarily reflect the results obtained in experimental approaches. In the latter, archaeal primer pair ArchV34 showed the highest similarity to the archaeal community structure compared to observed by the metagenomic approach and thus appears to be the appropriate for analyzing archaeal communities in biogas reactors. However, a disadvantage of this primer pair was its low specificity for the archaeal domain in the experimental application leading to high amounts of bacterial sequences within the dataset. Overall our results indicate a rather limited comparability between community structures investigated and determined using different primer pairs as well as between metagenome and 16S rRNA gene amplicon based community structure analysis. This finding, previously shown for Bacteria, was as well observed for the archaeal domain.